12 Oct 1997 



IFUP-TH 48/97 



hep-ph/9710312 



Two-loop QCD corrections to 
charged-Higgs-mediated b —> 57 decay 

Paolo Ciafaloni 

IFAE - Grup de Fisica Tedrica, Edifici Cn, Universitat Autdnoma 
de Barcelona, 08193 Bellaterra, Espana and INFN - Frascati - Italy, 

Andrea Romanino and Alessandro Strumia 

Dipartimento di Fisica, Universitd di Pisa and 
INFN, sezione di Pisa, 1-56126 Pisa, Italia 



Abstract 

The charged-Higgs-mediated contribution to the Wilson coefficient of the b — > 37 
magnetic penguin is expected to be one of the more promising candidates for a super- 
symmetric effect in B physics, probably the only one in gauge-mediated models. We 
compute the two-loop QCD correction to it. With naive dimensional regularization 
and MS subtraction, for reasonable values of the charged Higgs mass and for Ji = m t , 
we find a (10 -v- 20)% reduction of the corresponding one-loop effect. 

1 Introduction 

In supersymmetric models with sparticle masses mediated by supergravity Jl}, physics at the unification scale 
leaves its imprint in the supersymmetry breaking terms. As a consequence we expect that the sfermion mass 
matrices contain new sources of flavour and CP violation, cither due to unification physics [|| or related to 
the generation of the flavour structure itself [||. In this scenario it is quite possible that the effects due to 
virtual sparticle exchanges (mostly fj, — > ej, [i — > e conversion in atoms, electric dipoles of the electron and of 
the neutron, CP-violation in the -ftT-system and -B-systems), will be discovered even before than the sparticlcs 
themselves. 

Alternatively, in supersymmetric models, like Gauge-Mediated (GM) models, where the supersymmctry- 
breaking soft terms are generated at a lower scale where non-MSSM physics has decoupled, we expect that the 
only flavour violation present at low energy be the supersymmetrized extension of the standard CKM matrix. 
In this scenario the new supersymmetric flavour violating interactions (mainly the ones at charged Higgs and 
higgsino vertices, present in any realistic supersymmetric extension of the SM) are essentially unrelated to 
supersymmetry breaking and their flavour structure is controlled by the same CKM matrix. Consequently they 
do not introduce any new CP-violating phase nor they affect leptons, so that supersymmetry gives contributions 
only to 'standard' flavour and/or CP violating effects. These effects can nevertheless lead to detectable deviations 
from the expectations of the standard model in a few crucial observables in flavour physics, mainly in the b —> S7 
and b — ► s£ + £~ decays j|, |J. Supersymmetric effects in other processes are less interesting, because the SM 
background is larger and/or plagued by larger QCD uncertainties. 

From the point of view of the effective theory below the Fermi scale, in this scenario supersymmetry is 
expected to give a detectable correction to the Wilson coefficient, C7, of the b — s", magnetic penguin operator, 
usually named O7. This coefficient can be extracted from B.R.(6 — > S7) with a ~ 10% theoretical uncertainty || 
(jj. Alternatively it can be deduced from the spectrum of b — > s£ + £~ decays away from resonances. These decays 
have a lower branching ratio (around 6 • 10~ 6 for I = e, fi and smaller for £ = r) than b — » sj. Consequently, 
for a precise determination of CV(/2 rs m t ) one needs to wait for sufficient statistics. However, with data from 
5 • 10 8 BB pairs, the ler error on C7 will again be slightly larger than ±0.01 j§|. 
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In the SM the perturbativc QCD uncertainties in B.R.(6 — > sj) have been reduced from ~ 30% to ~ 5% 
computing the full Next-to-Leading-Order (NLO) QCD corrections. This computation can be divided in three 
steps: (i) at some scale /2 close to the electroweak scale, the couplings of the effective theory (containing only 
the light degrees of freedom) are determined up to 0(az) (|, in such a way that the full and effective theories 
describe equivalent physics; (ii) the effective theory is evolved via RGE techniques at NLO from the electroweak 
scale down to the B scale ||; (iii) at the B scale the matrix elements for the b — ► S7 process are computed 
with 0(0.3) precision With MS subtraction, and in the 'naive dimensional rcgularization' (NDR) scheme 
commonly employed, these three parts are of comparable numerical importance fl0|| . 

In the scenario where supersymmetry affects the low energy theory only giving extra contributions to the 
same Wilson coefficients that describe the SM physics, 'only' the step (i) is missing for a complete computation 
of supersymmetric corrections at NLO precision. We do this computation in the case of the charged Higgs 
mediated correction to C7. In the next section we explain why this correction is of particular interest. In 
section 3 we collect the ingredients for the computation. The computation in the full theory is presented in 
section 4, and the one in the effective theory in section 5. The final result is given in section 6. For a reasonable 
Higgs mass, ran ~ 2m t , we find that in the NDR scheme the one-loop charged-Higgs contribution, evaluated 
at p, = m t , is reduced by two-loop QCD corrections by 15%. 



2 Relevant supersymmetric effects 

Once that C7 will be measured with sufficient precision, the accurate SM computations could say that some new 
physics is required. Needless to say, however, as always in the case of radiative corrections effects, the discovery 
of non SM effect would not allow an immediate identification of its physical origin. The supersymmetric effects, 
in particular, depend on many unknown parameters^]: the charged Higgs contribution depend on m# and tan /3; 
the chargino contribution (at one loop) depend on the /x-parameter, on the weak gaugino masses, on the trilinear 
term of stop squark, A t , and on the up-type squark masses. So, it docs not seem useful to replace an approximate 
function of many unknown parameters with a more precise function (the NLO chargino contribution depends 
also on the masses of various supersymmetric particles in a non-decoupling way, if the couplings A t , At>, (72 at 
chargino vertices are expressed in terms of measured SM parameters). 

To discuss this point and to assess the relative importance between the different contributions we need to 
look in more detail at the one- loop predictions. In the simplifying limit uih S> m t , mg 3> m^ » ji > A-% 
the relevant contributions from Standard Model (SM), charged Higgs (H) and charginos (\) to the Wilson 
coefficients C7 of the b — > ,57 magnetic penguin, and of the b — > sg chromo-magnetic penguin are 

Cl M « -0.2 C 8 SM « -0.1 



2M% \36 tan 2 /3 3 m 2 2 J 8 ~ 2M% \Q tan 2 /3 m 2 2 J U) 



r x^ 1 (2 /i(A t tan/? + M ) mg 13 \ _ m 2 1 I 1 fx(A t tan/3 + /j) \ 

7 ~ 2m? sin 2 /3 ^9 + m? /i 2 6 ') 8 ~ 2m? sin 2 (3 \ 12 2 m 2 ) 

We see that the charged Higgs gives an important correction, with the same sign of the SM contribution. For 
example, even a heavy charged Higgs with mass run = 700 GeV gives a 10% enhancement of B.R.(& — > 57) 
over the SM prediction. The chargino contribution is relevant if a stop state is relatively light and the other 
up-squarks are heavier. A possible gluino/squark contribution is at most few % than the SM result, even in 
presence of new CKM-like mixing matrices with elements of the order of the CKM-oncs §, §. A possible 
neutralino/squark contribution is even more negligible. 

Moreover, the scenario where b — > 57 seems the most interesting candidate for a detectable supersymmetric 
effect can naturally be realized in gauge- mediation models [l2| . In these quite predictive models the chargino-up- 
squark loops do not give large contributions, unless tan /3 is large |l3| ] . This happens because in a typical gauge 
mediated spectrum the squarks are rather heavy (their mass terms are mediated by strong gauge interactions) 
and 'rather' degenerate among them (gauge interactions are generation universal), so that chargino/squark 
contributions are small and GIM-supprcsscd. Consequently, in gauge-mediation scenarios, the most (and, 
probably, only) interesting supersymmetric effect in B physics is given by a charged-Higgs / top correction to the 
b — > S7 magnetic penguin. This extra contribution depends only on mn and tan /3, and more precisely: 



^^We use standard notations for the supersymmetric parameters; in particular ran is the charged Higgs mass, tan/3 is the ratio 
between the vacuum expectation values of the two Higgs fields, and fi is the '/^-parameter', not to be confused with the MS scale p,. 
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• The dependence on tan/3 of the charged-Higgs- mediated magnetic penguin is very weak. As clear from ([I]), 
the term suppressed by l/tan 2 /3 is very small already for tan/?>2. We remember that smaller values 
of tan /?, close to one, are unnatural in the successful scenario of radiative breaking of the electroweak 
symmetry ]i4j , and make difficult to accommodate the large value of the top mass. 

• The charged Higgs mass run is not predicted by gauge-mediation scenarios; however the mass parameters 
of the Higgs potential determine the electroweak scale, so that (unless tan/3 is large), barring unnatural 
accidental cancellations, the charged Higgs should not be much heavier than the Z boson. 

So, the charged-Higgs mediated contribution depends almost uniquely on the charged Higgs mass and has 
defined sign. For these reasons we think that it is interesting to have a more accurate determination of the term 
not suppressed by 1/ tan 2 (3, that arise from graphs that contain one vertex A& b^t l H_ . 

More generally this computation also applies to two-Higgs doublet models, where the charged Higgs mediated 
contribution can be the only correction to the SM prediction. 



3 Preparation for the computation 

To be more precise, and to conform to the standard notations, we recall the effective Hamiltonian for the b — > S7 
decay 



i=l 



We have approximated V c bV* s a —VtbV t * s and, in the limit m s — > 0, 

O x = 4(s l ^V L c>)(c ] lti r h b l ) 2 ee l^V^cP^ViV) 

3 = 4(5*7^6*) E(g J 7M?W) 04 = 4(s V^') X^'VW) 

7 ee 4-^m 6 (s 7 ^R&)i^ 8 = ij^m b (s^T a P K b)G a ^ 

where i,j are colour indexes, and the sum is over all light quarks q = {b,c, s,u,d}. The relevant Wilson 
coefficients at the electroweak scale, in LO approximation, are 

0| M = 1 Cf = 

C 7 = 2^ uP FE + OwPbe] C 7 = -[{q u PFI + qHPBl) H tan2 -q ] (4) 

r SM _ 3 p /^H _ 1 [n 1 ^P-B 1 



2 2 L tan 2 /3 J 

At leading order the magnetic penguin coefficients CV and C% are given by the value of the corresponding 
penguin diagrams, without corrections from the matching procedure. The magnetic penguin loop-functions 
PF/B.i/E{ r )-, that we will also employ to write the divergent part of our two-loop Feynman graphs, are listed 
in the appendix. The indices F or B indicate if the external photon (or gluon) is attached to a Fcrmion or 
a -Boson. The indices E or / indicate if the hclicity flip factor, m&, comes from the -External b leg, or from 
an internal vertex of the graph. The parameter r = rn 2 B /m 2 F is the ratio between the (squared) masses of the 
boson B and of the fermion F in the loop (r = M^, / m 2 in the SM case and r = m 2 H / m 2 in the charged Higgs 
case). 

We will not compute the NLO corrections to the b — > sg chromo-magnetic penguin operator 0g . This would 
be formally inconsistent in an expansion where no small parameter is present in the QCD RGE corrections so 
that the RGE mixing between O7 and 0% is considered of 0(1). However the relevant RGE- loop factor is of 
order I = (a.3 /An) In M^/m 2 f» 0.1, smaller than 1. More in detail, at leading order, the coefficient C7 at the 
B scale is obtained in terms of the coefficients at the electroweak scale as 

C 7 (m b ) a U 72 + U 77 C 7 + U 7S C S a -0.155 + 0.7(-0.2 + Cf USY ) + 0.085(-0.1 + Cf USY ) a -0.30 + 0.7Cf USY 

Here Uij is the evolution matrix from the weak scale to the B scale for the Wilson coefficients CV This confirms 
that it is not necessary to compute 08 at NLO. 

More generally, the RGE-loop factors of order one, that surely need to be exactly resummed via the RGE 
techniques, are the ones that determine the running of the strong coupling constant. This happens because the 
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QCD (3 function receives contributions from all ('active') quarks, while only some specific flavours contribute, 
i.e., to the mixing between the b — > s penguins 7 and 0$. 

Consequently, when studying the decay b — > s£ + £~, it is not inconsistent to add the correction to C 7 we 
are going to compute (together with the corresponding low energy part), even if not required by a 'formal' 
expansion at NLO fll5| . 

This correction has already been computed in |]l6f with effective theory techniques, but only in the limit 
rriH 3> TUt where the effect is negligible. The complete NLO computation has instead been done in the case of 
charged Higgs corrections to the -Bo-Bo mixing fl7f . 

We can now pass to the NLO computation of the part of the charged-Higgs mediated contribution to C 7 not 
suppressed by 1/ tan 2 /3. This computation requires to match the NLO (two- loop) b —> amplitude in the full 
theory with the one of the effective theory. We choose to employ the equations of motion (idb = m^b) so that 
7 is the only relevant operator. In this way we are forced to match on-shell b — * S7 amplitudes, plagued by 
infrared divergences that must be properly treated. 

In both versions of the theory we choose the Feynman gauge for the gluon propagator. Since we never need 
to define traces like Ttj^j^ ■ ■ - "f p , nor the completely antisymmetric tensor, we can employ naive dimensional 
regularization (i.e. anticommuting 75) with MS renormalization scale p,. This regularization is more convenient 
than the supersymmetry-preserving 'dimensional reduction'. Infact the contribution we want to compute does 
not depend on supersymmetric parameters (like gaugino and sfcrmion masses,. . . ) that can be related among 
them making assumptions about physics at high energy, giving rise to predictions usually computed employing 
dimensional reduction. 

For simplicity in the following we denote with C 7l Cs,- ■ ■ , only their charged-Higgs mediated part that we 
want to compute. 



4 b — > S7 in the full theory 

As said we are interested in diagrams that contain one A& b^t^ H_ vertex. We have to compute the 16 two-loop 
diagrams shown in figure 1. 

We write the various contributions T to the two-loop b — > S7 effective Hamiltonian in the full theory as 

n b ^ = -^--V tb V t * s -(C ro + ^c 3 M)^-, M = ^M r (5) 

where C3 = 4/3 is the quadratic Casimir for the fundamental representation of SU(3) C , and the overall factor 
has been defined in such a way that the C3M is normalized as C71, the NLO term of the Wilson coefficient C7 
of the operator 7 in the effective theory: 

C 7 = C 7Q + ^C 71 + 0(aj). 

47T 



The LO coefficient of 7 in the full theory, C 7 q, differs from the corresponding one in the effective theory, O 



by terms that vanish as d —> 4, as described in eq. (A.l) in the appendix. Since the parameters of the theory 
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receive divergent corrections, these terms will give a finite contribution to C 7 \. 



4.1 Two loop diagrams 

A given Feynman diagram T gives a contribution Air = J Ft (where J denotes the standard two-loop integration 
over internal loop momenta) that depends on the heavy masses, M = {TOt,m#} and on the light masses and 
momenta, m = {to;,, m s ;pb,p s }, of the b and s quarks. The full expression Mr(M, to) contains negligible terms 
suppressed by powers of m/M. The natural trick that allows to avoid computing all the unnecessary terms, 
without loosing terms like In m/M, consists in adding and subtracting to the loop-integrand Fr(M,m) some 
appropriately chosen simpler term AFr(M, m) that has the same low energy structure of the full Fr(M, to). The 
result without irrelevant powers of m/M is given by the simpler integral J[Fr(M, 0) — AFr(M, 0)] + AFr(M, m). 
In practice simpler means that AFr can be chosen as a product of two one- loop integrals, instead of a two-loop 
integral, so that, in dimensional regularization, J AF(M, 0) = 0. 

It is convenient to employ a standard technique |lq| , usually named 'Heavy Mass Expansion' (HME), that 
generalises and systemizcs this procedure. It is possible to show that, in dimensional regularization, the ex- 
pansion of a Feynman diagram T, that depends on 'heavy' masses and momenta M and on 'light' masses and 
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Figure 1: The 16 two-loop Feynman graphs (we show on a single diagrams all possible attachment of the photon 
and write near to it the name we give to the diagram) that we have to compute. The symbol ® denotes the 
vertex bptLH- . 



momenta to, in the limit M 3> m can be done at the level of integrands as 

F r (M, m) ~ F r/ 7 (m) T Qi ^ F 1 {q 1 ,m 1 , M) asM»m (6) 
7 cr 

where T x indicates Taylor expansion in the variables x up to the desired order. The sum is performed over all 
the "asymptotically irreducible" subgraphs 7 of T , i.e. those which satisfy the following two conditions: 

1. 7 contains all the lines with heavy masses; 

2. 7 consists of connectivity components that are one-particle- irreducible with respect to lines with small 
masses. 

Finally q 1 denotes the set of momenta external to the subgraph 7. _F 7 is its Feynman integrand, and fr/7 is the 
Feynman integrand of the reduced graph r/7. We refer the reader to |]lo| , [lq| for explanations and examples. 
Here we discuss how this technique works in our case, and its relation with the natural trick previously discussed. 

- The case 7 = T is always present and gives a 'naive' Taylor expansion in m. This contribution corresponds 
to the term Fp(M, 0) of our previous example. 

- There is one non-vanishing contribution from the 'non naive' part of the HME expansion from each one 
of our graphs in fig. 1 (except the "GI" diagrams with an 'internal' gluon loop). The "asymptotically 
irreducible" subgraph 7 coincides with the Higgs/top loop that gives the LO result. This term of the HME 
expansion can be seen, at a diagrammatic level, as the contribution of the full diagram with the heavy 
propagators contracted to a point, and corresponds, in our previous example, to the term AFr(M, m). 
When computing it we neglect powers of p = m1/m%. 

- Terms analogous to — AFp(M, 0) are not included in this form of the HME, since they vanish in dimensional 
rcgularization. For this reason some fake infrared (IR) divergences can (and will) appear in the 'naive part' 
where m = 0. These divergences are cancelled diagram by diagram by fake ultraviolet (UV) divergences 
that appear in the non naive part of the HME expansion. 

The 'naive part' of the expansion also contains the usual UV divergences that cancel upon renormalization. 
Similarly, the 'non naive part' also contains IR divergences, cancelled in the final physical result after phase- 
space integration and inclusion of QCD bremsstrahlung b — ► sjg. 

In order to check in detail all these cancellations, we have separated IR from UV divergences. It has been 
possible to do it in dimensional rcgularization, since our graphs have only single poles 1/e. In the following we 
will denote as l/e uv an ultraviolet pole, and as l/£i r a pole of infrared origin. 
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Table 1 : Coefficients of the divergent parts of the single graphs. The loop functions Pbi and Pfi are defined in 
the appendix; the 'names' of the graphs are defined by fig. 1. 

In our case, since we have already a mj, factor from the vertex, the Taylor expansions have to be performed 
up to first order (up to second order in the 'non-naive part' of the diagrams 'FLIGE' and 'FRIGE' with four 
light propagators). The 'naive' parts of diagrams that differ by the exchange b # s (like the diagrams 'FLIGE' 
and 'FRIGE') are equal because the Dirac equation, that sees the difference between the b and s masses, is 
never used. 

We have written a Mathematica code that simplifies the spinor algebra, expands the 'light' factors in the 
loop integrands in Taylor expansion up to the appropriate order, reducing the full expression to a sum of one 
and two-loop scalar integrals that have been separately computed. This technique is less efficient than the 
alternative one based on Feynman parameters; however it is easy to teach it to a computer. The non-naive part 
of the HME expansion is computed in the same way, performing the Taylor expansion also in one (appropriately 
chosen) loop momentum. In this way the whole computation is done in few minutes by a normal 'personal' 
computer. We show the separate result of the two parts of the expansion (A^2 loop from the 'naive' (high-energy) 
part, and A^hme from the 'non naive' (low energy) term of the HME expansion) in eq.s (|l^). We have listed 
in table 1 all the divergences present in the single graphs. 

4.2 Renormalization 

It is necessary to rcnormalizc the top mass, m t , the bottom mass, nib, and to express the result in terms of 
canonically normalized b and s quark fields. 

The simplest way of renormalizing the top mass consists in expressing the bare top mass, mto, in term of 
the renormalizcd 'running' mass, m t , 

5m t a 3 3 

mt = m t0 + dm t , = c 3 - 

mt 47r e uv 

(cs = 4/3) in the one loop contribution, CVo(r). The corresponding contribution A4 t to the amplitude M is 

M t = —2rC' m (7a) 

where ' denotes derivation with respect to r = m 2 H / m\ . The pole 6-quark mass nib is expressed in terms of the 
bare b mass, mi,o, as 

Sm b a 3 3 fl 2 

m b =m b0 + dm b , = c 3^— ( h31n^- + 4). 

m b 47r e uv mi 

Since in the one-loop graphs the m b Q factor comes only from the vertex (we never use the Dirac equation to 
reduce the operator basis), the contribution Ai b to A4 from the renormalization of m b is simply given by 
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Figure 2: The Feynman graph that accounts for 
the renormalization of the mass and kinetic terms 
that mix the b and s quarks. 



3 n 2 

M b = -C 70 (- +3^-^+4). (7b) 



e m b 

The counterterms for on-shell renormalization of the wave-function of the b and s quarks are 

"3,1 , 2 , o p? 



Z q = l-c 3 -l( + — + 31n^ + 4), q = {s,b}. 

47T £uv Sir J7l„ 



Expressing the one-loop result in terms of canonically normalized b and s quark fields, we get the following 
correction Mz to M. 

M z = -C 70 (- + — + 3 In -^+4). (7c) 
e e-„ m b m s 

Finally, top/Higgs loops generate non diagonal kinetic and mass terms between the b and s quarks that is 
necessary to eliminate via appropriate flavour rotations of the quark fields. This gives rise to a contribution to 
M. when a magnetic bbj dipole operator is rotated into a bsj dipole. However, it is more convenient to take 
into account this particular correction adding the (one+one)-loop diagram of figure 2 (the similar graph with 
b ^ s is suppressed by powers of m s /m,b) to the list of the two-loop diagrams. 
We denote the renormalization contributions as M c tr = Mt + Mb + Mz- 



5 b — > s-f in the effective theory 

Since all the supersymmctric contributions (if i?-parity is conserved) arise first at one loop level, their mixing 
structure is simpler than the one of the SM contributions. In particular all the LO supersymmctric contribu- 
tions to the Wilson coefficients C, are given by the corresponding penguin diagrams without extra matching 
contributions. We list here all the various contributions to the b — > S7 decay amplitude in the effective theory 
relevant for our computation (we are neglecting terms suppressed by tan -2 /?). We denote by M the value of 
the on-shell b — > sj amplitude in the effective theory, normalized in the same way as M (the corresponding 
amplitude in the full theory). The non zero contributions to M are, in our case: 

• The contribution A4j given by the matrix element of O-j itself at order 0:3 



Mr = — C71 

C3 



C70 In p 



h m y 

£ir ml 



2 In p 

2 H 



(8) 



The first term is the tree- level matrix element of the term we ultimately wish to extract. The other terms 
are obtained by the one-loop QCD correction to the leading order magnetic penguin O7. 



• The order 0:3 matrix element of the chromo-magnetic penguin operator 0% 

-,2 



M s = q d Ci 



81] 



4 

^UV 



4 In 



11 



2tt2 
3 



2in 



(9) 



• As the operators mix under renormalization, we have to consider countcrterm contributions induced by 
operators of the form Ci SZij(sj\Oj\b). Using the known renormalization constants ]l9| ], the non- vanishing 
contributions to the b — > S7 amplitude M are 

4 4 
M77 = C 70 , and Ms7 = qd C &0 . 

Suv ^uv 



'Evanescent' operators do not give contributions. 
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• Like in the full theory, is of course necessary to renormalizc also the parameters and holds of the theory, 
that in this case are the 6-quark mass and the b and s wave-functions. The counterterm due to the 6-quark 
mass rcnormalization is 



M b = -O 



70 



_i_ +3m iL +4 



(10) 



when using the pole 6-quark mass as in the full theory, while the rcnormalization of the b and s fields 
gives a contribution Mz 



M z = -O 



70 



1 2 nl p 2 
— + — +3In-i 

e uv m b m s 



(11) 



The amplitude M in the effective theory is UV finite and has an IR divergence — (2 + h\ p)Cm/ 

6 Matching and final result 

We now collect all the necessary terms and illustrate the matching procedure. The amplitude in the effective 
theory is 



M = —C-i + q d C 80 

C3 



2rr 2 , p 2 
11 — 2i?r — 4 In At 



-G 



7(1 



(2 + Inp) {p2,mlr + 8 + 4 In 4 + l^p 1 In 2 p 
£i r mf 2 2 



The amplitude M in the full theory is the sum of the contribution Mi loop, defined as the 'naive' part of the 
two-loop diagrams, of the 'non-nai've' part, A^hme, and of the counterterms, M c tv 
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M ctl - 
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6rC 70 - 4C 70 



-G 



7(1 



2 (A7K)! + 8 + 41n 4_3 

£ir ml 2 



(12a) 
(12b) 
(12c) 



where the two-loop function f(r) is, for arbitrary electric charges q d = qu + q u 
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(13) 



We have neglected higher powers in e and p = m 2 /m 2 , and we have included in A^hme the diagram in figure 2. 
The full amplitude M is 



/(r) + (4C 70 - 6rC 70 + 4q d C S0 ) In A, + ^C, 



2tt 2 fl' 
— -Il-2i7r-41n — 
3 m 



-C70 



(2 + Inp) 



(M 2 /m 2 )^ 



M , 1 



1 



4 In — 2 H — m P l n /° 



m 6 2 



0(e). 



The on-shell 6 — » S7 amplitudes (.M in the full theory and yV( in the effective theory) are infrared-divergent. The 
correct matching is obtained requiring that the infrared safe decay rate obtained including QCD brcmsstrahlung 
b — > S7<7 be the same in both (full and effective) descriptions. This gives jlO|, pT[ 



where 



Mb = (2 + lnp)G 



70" 



M + Mir = M + M 
(p/m b ) 2 



M iT = (2 + lnp)C- 



70" 



{p/m b f 



(14) 
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So, the final result for the NLO charged Higgs contribution is (setting qjj = — 1 and q u = 2/3) 



C3 



/(r) + (4C 70 - 6rC 70 + 4q d C 80 ) In 



4 2 - 13r + 7r 2 2r 7 - 64r + 33r 2 



(r - l) 3 



9 



16 2 - 7r + 3r J 
~9 (r-1) 3 



Lia(l-r) 



2 8 - 47r + 21r 2 4r 3 + 14r - 8r 2 



(r - l) 2 



9 (r-1) 4 



-1)4 
In r I In — ^ 



In the limits m# = mt and to// 3> mt the charged Higgs contribution becomes 
C 7 (r-1) = 



CV(r — ► 00) 



36 4tt 
1 lnr 03 1 
4r 3r 47r 9r 



181 35 , p 2 
In 



81 27 " m 2 



2(42 + 4tt 2 - 33 lnr + 12 In r) + (42 - 32 In r) In ■ 



In ] 



(15a) 

(15b) 
(15c) 



For ma = 500 GcV and // = m t the term we have computed gives a —17% correction to the LO value. 



7 Conclusion 

In a class of supersymmetric scenarios, i.e. in gauge mediated models, the charged Higgs mediated contribution 
to the b — > 57 magnetic penguin is expected to be the most interesting supersymmetric effect in i?-physics. 
The NLO QCD correction to its Wilson coefficient CV(/z), that we have computed employing naive dimensional 
regularization with MS subtraction, gives a —(10 -j- 20)% correction to the corresponding LO result (for p = m t 
and realistic charged Higgs masses). 



Note added: Together with our article, another one containing the same computation 20 has recently 
appeared. Our results are in perfect agreement, even if the conclusions are apparently different. In particular 
the authors of claim that, in the usual two-Higgs doublet model, the inclusion of NLO effects enhances the 
B.R.(& — ► S7) with respect to its LO value. This is due to the already known NLO enhancement in the SM 
contribution H [?], [l^]. On the contrary, the new correction to CV(/2 w m t ) that we have computed tends to 
reduce the charged-Higgs mediated correction. As discussed in the text, we consider CV(/I s=s mt) a more useful 
phcnomenological quantity than B.R.(6 — > S7). 
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A Useful functions 

The penguin one-loop functions P(r) employed in cq.s (Q) are given by P — P| £ ^o- where 

, p, 2 .r 2 - 1 - 2rlnr 3(r 2 - 1) - 2r(2 + r) lnr + 2rln 2 r 
Pbi = (1 + eln— )- 



Pfi = (l + eln^o) 



2(r-l) 3 4(r-l) 3 
p? 1 - 4 r + 3 r 2 - 2 r 2 In r 1 - 8 r + 7 r 2 - 6 r 2 In r + 2r 2 In 2 r 



m 



V 2(r-l)3 4 (r-1) 3 



, p 2 2 + 3r-6r 2 +r 3 + 6rlnr 22 + 27r - 54r 2 + 5r 3 - 6r(-6 - 6r + r 2 ) lnr - 18rln 2 r 
Pbe = 1 + eln-^ — — +e — ^ — — 4 '- 

mf 12(r - l) 4 72(r-l) 4 

, p 2 , 6r - 1 - 3r 2 - 2r 3 + 6r 2 In r 5 - 54r + 27r 2 + 22r 3 - 6r 2 (9 + 2r) In r + 18r 2 In 2 r 
Pfe = (l + eln-^o) 



12(r-l) 4 72(r-l) 4 

The charged-Higgs mediated one-loop result for the coefficients C70 of the b — » sj magnetic penguin, and Cso 
of the b —> sg chromo-magnetic penguin are (J, ^ 

C70 - \{q H P B i + quPpi) + 0(e 2 ), Cso = \ppi + 0(e 2 ). (A.l) 
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The bi-logarithmic function, hi2(x), employed in eq. ( fi"3| ) is defined as 

Jo s 
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